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ABSTRACT 

The  thickness  of  weak  shock  fronts,  across  which  the  pressure 
difference  was  in  the  vicinity  of  ten  atmospheres,  was  measured  by 
the  reflectivity  method  in  acetone,  ethyl  alcohol  and  ethyl  ether* 
They  are  all  of  the  order  of  10“^  cm,  thick.  Lower  limits  to  the 
thickness  were  established  in  water,  carbon  tetrachloride,  benzene 
and  ethylene  glycol,  A theoretical  expression  for  the  thickness 
is  derived.  Comparison  between  theory  and  experiment  suggests  that 
the  effective  viscosity  in  the  shock  front  is  somewhat  greater  than 
the  ordinary  shear  viscosity. 
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It  has  proved  possible  to  investigate  the  transition  region  in 
shock  fronts  in  gases  by  making  use  of  the  fact  that  the  thickness 
of  the  transition  region  is  of  the  same  order  of  magnitude  as  the 
wave-length  of  visible  light  under  suitable  conditions  of  shock 
strength  and  initial  gas  pressure;  because  of  this  circumstance 
the  optical  reflectivity  is  sensitive  to  the  thickness  and  density 
profile  of  the  shock  front. 

1,  D.  F.  Hornig,  Phys.  Rev.  179  (1947) 

From  their  reflectivity  the  thick- 
ness of  shock  fronts  has  been  measured  in  argon,  nitrogen,  oxygen 
and  hydrogen  and  rotational  heat  capacity  lags  have  been  studied  in 
the  latter  three  molecules. 

2,  G.  R.  Cowan  and  D.  F,  Hornig,  J.  Chem.  Phys.  18,  1008 
(1950) 

3.  E.  F.  Greene,  G.  R.  Cowan  and  D.  F,  Hornig,  J.  Chem.  Phys. 
19,  427  (1951) 

4.  E,  F.  Greene  and  D.  F.  Hornig,  J.  Chem.  Phys.  617  (1953) 
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Very  little  has  been  published 

(5) 

relating  to  the  thickness  of  shock  fronts  in  liquids.  Cole 

5.  R.  H.  Cole,  ’’Underwater  Explosions”  (Princeton  University 
Press)  Princeton,  N.  J.  1948,  p.46. 

has  published  a summary  of  the  work  in  underwater  explosions,  and 
in  it  he  states:  ”A  number  of  writers  have  made  such  estimates  (of 
the  thickness  of  intense  shock  fronts,  based  on  hydrodynamic  theory) 
....  the  common  result  is  ....  that  they  predict  thicknesses  of  the 
order  of  10”^  or  10“  cm  No  detailed  theoretical  work  has 

been  done  to  date  on  this  subject,  and  there  has  been  no  experimental 
evidence  as  to  the  thickness  of  any  shock  in  a liquid.  In  the 
present  paper,  the  thicknesses  of  weak  shocks  in  liquids  will  be 
calculated  from  the  hydrodynamic  equations  and  an  equation  of  state. 
An  apparatus  to  measure  these  thicknesses  will  be  described,  and 
experimental  results  in  a variety  of  liquids  will  be  presented. 

Shock  Front  Thickness  in  Liquids. 

The  three  hydrodynamic  equations  representing  the  conservation 
of  mass,  momentum  and  energy  across  a plane  shock  front  progressing 
in  the  positive  X direction  in  a homogeneous  compressible  fluid 

may  be  written  respectively 

jOU.  - a,  (1) 

j /^5jT  = ^ (2) 

K ^ = K.  E ■<- fc> — e 

where  p is  the  density  at  any  point  in  the  fluid,  u the  particle 
velocity  with  respect  to  the  shock  front, ^ the  shear  viscosity, 
p the  pressure,  E the  internal  energy  per  unit  volume,  K the  coeffi- 
cient of  thermal  conductivity  and  T thoaabsolute  temperature, 


3. 


B-t  t and  c are  constants  defined  by  the  initial  and  final  condi- 
tions. A convenient  equation  of  state  is  obtained  by  expanding  the 
pressure  as  a function  of  the  density  and  the  entropy: 


Weak  shook  fronts  in  liquids^across  which  the  pressure  differ- 


ence Ap  is  of  the  order  of  ten  atmospheres  are  to  a good  approxi- 
mation isentropic.  Therefore,  the  terms  in^S  in  Eq.  (4)  may  be 
dropped  and  only  those  in  and  <^^^used.  This  approximation  is 
equivalent  to  replacing  Eq,  (3)  by 


(5) 


The  thickness  of  the  shock  front  is  defined  by  the  expression 

ex- Pi 


(6) 


YY>e>.x  ■ 


where  the  subscripts  1 and  2 refer  to  the  limiting  regions  ahead 
of  and  behind  the  shock  front,  Eqs,  (1),  (S)  and  (4)  may  be  solved 
for  dp/dx  to  yield,  in  the  case  of  weak  shocks. 

/ — T 


In  deriving  this  equation  use  has  been  made  of  the  fact  that 
aui  + pi  - b = 0 in  the  limiting  region  before  the  front,  and  that 
U]^,  the  velocity  of  the  shock  front,  differs  from  the  sound  velocity 
Cq  by  less  than  one  part  in  one  hundred.  Use  was  also  made  of  the 
relations 
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and 


(8) 


^ W 

' P/  ^ P 

The  last  two  terms  In  in  Eq.  (7)  are  actually  third  order 

as  a consequence  of  Eq,  (9) , For  example,  for  a 10  atmosphere 

shock  in  ethyl  alcohol,  they  amount  to  less  than  one  part  in  700 

of  the  complete  term  in  if  the  coefficients  are  taken  from 

(6)  * 

the  data  of  Amagat, 


6.  E.  H.  Amagat,  Annales  de  Chimie  et  Physique  505  (1893) 

Consequently,  to  second  order  eqn.  (7) 


becomes 

-Ji  — ifl 

3 f cu 


so  that 


h 


L = 


C ^ 
_ 

(10) 


The  coefficient  has  not  been  determined  for  any  of 

the  substances  studied.  It  could  best  be  obtained  from  a study  of 
the  variation  of  the  velocity  of  sound  with  pressure,  and  tempera- 
ture, or  with  density  and  temperature,  since  it  is  given  by  either 
of  the  following  two  relations: 


(11) 


( ijl  J ,5  Co  ( 1^1  + L___.  / , Co\ 
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where  In  addition  to  the  symbols  previously  defined,  oi.  is  the 
coefficient  of  thermal  expansion  and  is  the  compressibility. 
Sufficient  data  to  use  these  relations  were  not  available  in  the 
literature. 

However,  as  a rough  approximation  to  one  may  use 

Cp/Cv  available  from  the  compressibility 

measurements  of  Amagat.  In  doing  this  it  should  be  noted  that 
Cp/Cv  taken  from  this  source  differs  from  by  11^  for 

ethyl  alcohol,  20%  for  acetone  and  36%  for  ethyl  ether.  There  may 
also  be  corresponding  errors  in  the  second  derivative.  The  ex- 
pected shock  front  thickness  is  estimated  for  a number  of  substances 
in  Table  I. 
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Theoretical  Shock  Thickness 
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In  Interpreting  Table  I It  must  be  realized  that  in  using  the  ordi- 
nary shear  viscosities  no  account  has  been  taken  of  the  possibility 
that  ytA  depends  on  the  rate  of  shear#  Because  compression  in  the 
shock  is  equivalent  to  sonic  processes  of  the  order  of  300  mega- 
cycles/sec,, this  is  not  necessarily  a safe  assumption.  All  of  the 
values  listed  in  Table  I were  taken  at  20°C,  except  the  compress- 
ibility data  for  acetone.  This  was  available  only  at  0®C, 

The  Effect  of  Bulk  Viscosity  on  the  Thickness 

In  recent  work  in  the  field  of  ultrasonic  absorption,  the 

validity  of  using  only  the  ^ear  viscosity  y>(_  to  describe  the  dissi- 

9 10  11  12 

pative  forces  in  a dynamic  system  has  been  questioned,  ’ * * 


9.  L.  Tisza,  Phys,  Rev.  ^1,  ^31  (1948) 

10.  Carl  Eckart,  Phys.  Rev.  7^,  68  (1948) 

11,  J,  J.  Markham,  R,  T,  Bey«r  and  R,  B,  Lindsay,  Rev, 
Mod,  Phys.  353  (1951) 

12,  L.  Hall,  Phys.  Rev,  73,  775  (1948) 


Tisza  has  proposed  that  in  addition  to  ^ , a bulk  viscosity  may  be 
effective  undor  dynamic  conditions.  If  k is  the  bulk  viscosity, 

Eq,  (11)  becomes 


4 Co 


/ill 


L c3c/V 

. V^-  J 

(12) 


If  the  value  of  k/^  is  of  the  order  of  3,  as  has  been  roported^^ 
for  some  liquids,  L*  would  be  aboTit;  3L,  Accurate  experimental 


measurements  of  L combined  with  a well  established  value  of 
should  yield  informfition  as  to  the  bulk  viscosity  forces  acting 
in  the  shock  process. 
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Experimental  Method 

In  order  to  apply  the  reflectivity  method,  a knowledge  of  the 
variation  of  index  of  refraction,  n,  with  distance  through  the 
shock  front  is  necessary.  This  may  be  obtained  by  integrating 
Eq,  (10)  to  obtain 

f,  ^ 


_ H 


Vl 


(13) 


} -h  e.  ■<--  j 

and  since  is  a linear  function  of 

= t-  . ciii.) 

It  is  interesting  to  note  that'  this  is  exactly  the  expression 
used  by  Cowan  and  Hornig  as  their  "Model  II"  function.  The  re- 
flectivity equation  is  therefore  the  same  as  that  of  Cowan  and 


Hornig: 


T^Lco5  9 
'h 


. f a. 


IT  COS  9 


(15) 


In  the  above  equation  R is  the  reflectivity,  ^ the  wave-length 
of  the  incident  light,  & the  angle  of  incidence  and  n-j_  the  index  of 
refraction  of  the  medium,  n is  the  total  change  in  index  of  re- 
fraction across  the  shock  front.  The  part  of  the  equation  enclosed 
in  square  brackets  is  plotted  as  a function  of  X/Z.co5^  in  Pig,  1, 

Two  essentially  different  methods  of  determining  L from  Eq,  (l5) 
suggest  themselves.  First,  following  Covran  and  Hornig,  the  parame- 
ter Loos  may  be  varied  bj  changing  either  X or  0 , and  the 

slope  of  the  reflectivity  function  determined.  The  value  of  L may 
then  be  detei’mined  from  Pig,  1,  Alternately,  the  apparatus  may  be 
calibrated  in  terms  of  absolute  reflectivity  and  Eq,  (l5)  solved  for 
L,  In  this  case  h is  given  by  the  Lorontz-Lorenz  formula  which 


is,  for  sufficiently  weak  shocks. 


9, 


A,  photograph  of  the  apparatus  used  is  shown  In  Pig.  2,  and  a 
schematic  diagram  in  Pig,  3*  The  apparatus  was  a modification  of 
that  used  by  Cowan  and  Hornigc  A shock  wave  was  generated  in  ni- 
trogen by  means  of  a bursting  diaphragm  in  a vertically  mounted, 
2-7/8”  diameter  shock  tube*  Copper  diaphragms,  ,008"  thick,  were 
used  to  produce  shocks  with  a pressure  difference  of  10  atmospheres. 
The  shock  wave  travelled  down  the  tube  until  it  impinged  at  normal 
incidence  on  the  surface  of  the  liquid  in  the  bottom  of  the  tube, 

A ,001”  membrane  of  duPont  ’’Mylar”  film  was  located  about  a centi- 
meter below  the  surface  of  the  liquid  and  served  to  separate  con- 
taminated liquid  above  from  the  clean,  carefully  filtered  sample  in 
a cell  at  the  base  of  the  tube.  The  3 cm,  layer  of  liquid  above  the 
cell  tended  to  slow  down  copper  fragments  and  protect  the  membrane, 
Windows  in  opposite  sides  of  the  cell  allowed  the  collimated  light 
beam  from  a carbon  arc  to  enter  the  liquid  at  a known  angle  to  the 
vertical;  a type  931A  photomultiplier  in  conjunction  with  an  optical 
system  received  the  reflected  light  pulse  as  the  shock  travelled  past 
the  window.  Suspended  matter  in  the  sample  had  to  be  carefully 
eliminated  since  it  caused  intolerable  light  noise;  in  fact,  the 
reduction  of  scattered  light  was  one  of  the  major  experimental  prob- 
lems in  this  work.  Samples  were  filtered  through  three  fine  sin- 
tered glass  filters  directly  into  the  cell,  which  was  constructed 
entirely  of  stainless  steel,  to  eliminate  suspended  particles  as  far 
as  possible. 


10, 


The  shock  tube  was  adjusted  to  within  2 minutes  of  arc  of  the 
vertical  by  means  of  a plumb  bob,  and  the  optical  zero  was  fixed 
within  30  seconds  of  horizontal  using  surveying  instruments.  Accu- 
rate angle  settings  were  imperative  since  the  angles  of  incidence 
used  were  close  to  grazing,  where  the  (1  + tan^^  ) factor  in  Eq.(lO) 
changes  rapidly  with  angle,  ^hock  pressures  were  computed  from  the 
measured  pressures  in  the  compression  and  expansion  chambers,  using 
previously  published  relations. 


13*  P.  W,  Geiger  and  0.  W,  Mautz,  ’’The  Shock  Tube  as  an 

Instrument  for  the  Investigation  of  Transonic  and  Super- 
sonic Plow  Patterns”,  U,  S.  Navy  Department  of  Naval 
Research  Contract  No.  N6-0NR-232  Task  Order  IV,  Engi- 
neering Research  Institute,  University  of  Michigan, 

June,  1949,  p.l2. 


A correction  was  made  for  the 

decrease  in  shock  strength,  due  to  non-ideal  bursting  conditions, 
reported  by  Bleakney  et  al,^^ 

14.  W.  Bleakney,  D,  K,  Weimer  and  C.  H.  Fletcher,  Rev,  Sci. 
Instr,  007  (1949). 


increase 


15 


as  well  as  for  the  more  than  two-fold 


l5.  R.  Courant  and  K,  0.  Friedrichs,  “Supersonic  Plow  and 
Shock  Waves”,  (Interscience  Publishers,  Inc,,  New  York, 

1948)  P.  1^3. 


in  shock  pressure  caused  by  reflection  at  the  liquid-gas 

Interface, 

An  image  of  the  carbon  arc  light  source  was  cast  on  a horizontal 
slit,  which  then  became  the  effective  light  source  for  the  rest  of 
the  system.  The  slit  was  designed  so  that  the  divergence  of  the 
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beam  in  the  vertical  plane  was  limited  to  plus  or  minus  30  minutes 
of  arc.  The  optical  system  was  arranged  so  that  the  direct  light 
beam  passed  through  the  cell  and  left  through  the  exit  window;  it 
was  absorbed  on  a black  paper  mask  outside  the  tube*  This  arrange- 
ment was  found  to  materially  reduce  the  background  of  scattered 
light,  A lens  and  mask  in  the  exit  system  permitted  only  nearly 
parallel  light  to  reach  the  photomultiplier,  still  further  reducing 
the  unwanted  scattered  rays.  The  entire  optical  system  was  mounted 
on  a 500  lb,  concrete  block,  entirely  separate  from  the  shock  tube, 
to  prevent  vibrations  from  reaching  the  detection  system. 

The  scattering  level  in  the  cell  was  measured  before  each  shot 
by  setting  the  angle  of  incidence  at  a standard  value,  and  measuring 
the  signal  obtained  when  the  arc  beam  was  ’’chopped”  by  a rotating 
sector.  Whenever  the  scattering  level  increased  following  a shot  a 
pin  hole  was  found  in  the  membrane.  High  scattering  levels  gave  rise 
not  only  to  random  ’’light”  noise  but  also  to  well  defined  spurious 
signals. 

Effective  angles  of  incidence  in  the  liquid  were  computed  by 
using  the  knovm  divergence  of  the  light  beam  in  connection  with  the 
reflectivity  equation  to  obtain  a properly  weighted  average,  after 
allowance  had  been  made  for  I’efraction  at  the  windows.  The  external 
angles  of  incidence  used  were  80®  and  82®,  corresponding  to  actual 
effective  values,  in  acetone  and  ethyl  alcohol,, of  82®45’  and  84®l5' • 
In  ethyl  ether  the  effective  values  were  82°37'  and  84®05’. 

Wavelengths  were  isolated  at  two  points  in  the  visible  spectriom 
by  Wratten  filters;  effective  wavelength  values  were  fixed  by  com- 
bining filter  transmissions,  measured  using  a carbon  arc  source  in 
connection  with  a recording  spectrograph,  with  the  theoretical 
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reflectivity  curve  to  obtain  a weighted  average  value.  The  detector 

In  the  spectrograph  was  a type  1P21  photoinultiplier,  which  has  th© 

same  spectral  sensitivity  as  the  93I-A*  Wratten  filter  No*  33' 

o 

yielded  a weighted  average  wavelength  of  4300A,  based  on  a shock 
thickness  of  10''^cm.  Wratten  No.  21  yielded  an  average  value  of 
570oX«  These  two  filters,  as  well  as  the  white  light  average  of 
473'oS,  were  used  in  all  the  experimental  work. 

Output  signals  were  displayed  on  a cathode  ray  oscilloscope  and 
recorded  photographically.  An  expanded  sweep  allowed  accurate  tim- 
ing measurements  to  be  made.  The  frequency  response  of  the  elec- 
tronic system  was  constant  to  100  kc.,  falling  off  to  $0%  at  400  kc. 

A signal  was  observed  in  some  liquids  upon  passage  of  the  pri- 
mary shock  wave,  A second  signal  was  observed  as  the  shock  was  re- 
flected from  the  bottom  plate  of  the  cell  and  again  traversed  the 
light  beam.  These  signals  were  identified  as  true  reflections  by 
their  timing,  their  shape,  their  duration  and  especially  by  their 
wavelength  dependence.  The  signals  were  observed  to  be  stronger  at 
longer  wavelengths,  eliminating  any  possibility  that  they  might  be 
caused  by  any  spurious  scattering  phenomenon. 

The  recording  system  was  calibrated  in  terms  of  absolute  re- 
flectivity by  using  a shock  front  of  known  thickness  in  nitrogen  as 
a reference,  A shock  wave  which  had  been  thoroughly  studied  by 
Greene  et  al^  was  generated  in  the  tube  by  duplicating  their  condi- 
tions, The  thickness  of  this  shock  was  known  to  be  1,63  x 10"^  cm, 
plus  or  minus  its  reflectivity  was  established  through  the  use 

of  the  Rankine-Hugeniot  relation  in  conjunction  with  Eq,  (l3)  and 
(16),  X /ir^LcoJ>  9 for  this  shock  falls  to  the  right  on  Pig,  1, 
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where  the  reflectivity  varies  slowly  with  thickness,  thus  enhancing 
the  accuracy  of  the  calibration.  The  thicknesses  observed  in  liquids^ 
on  the  other  hand,  correspond  to  points  on  the  left  of  Pig,  1,  in  a 
region  of  rapidly  varying  reflectivity,  which  tends  to  increase  the 
reliability  of  the  results. 

The  scatter  of  reflectivity  values  obtained  from  gas  calibra- 
tion shocks  was  considerably  less  than  that  observed  in  the  liquid; 
this  was  originally  attributed  to  the  complications  introduced  by 
the  passage  of  the  shock  through  the  liquid-gas  interface  as  well  as 
the  factors  mentioned  in  the  preceeding  paragraph.  However,  the 
condition  of  this  interface  at  the  instant  of  arrival  of  the  shock 
wave  was  investigated  by  totally  reflecting  a light  beam  from  the 
surface  and  observing  any  deflections.  No  disturbances  were  ob- 
served in  the  central  portion  of  the  surface  prior  to  the  arrival  of 
the  shock,  although  there  was  some  evidence  of  deformations  near  the 
edges. 

Experimental  Results  , 

The  experimental  results  at  26^3°^  are  listed  in  Table  II, 

Shock  wave  thicknesses  were  successfully  measured  in  acetone,  ethyl 
alcohol  and  ethyl  ether,  but  no  reflections  could  be  observed  in 
water*  carbon  tetrachloride,  benzene  or  ethylene  glycol,  even  from 
the  strongest  shock  waves  ( ^p  = 16,7  atmospheres)  available.  The 
minimum  reflectivity  observable  with  this  apparatus  was  about  RslO”"^; 
this  fefct  establishes  a lower  limit  on  the  shock  front  thickness  in 
these  liquids,  as  indicated  in  the  table.  The  thickness  values 
listed  in  Table  II  represent  the  average  of  1|  or  5 reflectivity 
measurements  at  each  of  the  angle  or  wavelength  settings,  except 
as  noted.  Probable  limits  were  computed  from  the  scatter  in  experi- 
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TABLE  II 


Experimental  Shook  Thicknesses  at  26“C* 


Liquid 

Acetone 


Shock  Strength 
(atmospheres) 


Ethyl  Ether 


Water 


Benzene 


13,1 


Ethyl  Alcohol  13.1 


16.7 


16.7 


Carbon 

Tetrachloride  16.7 


16.7 


Method 

Angle  dep. 
Angle  dep. 
Absolute  mag. 
Absolute  mag. 
Wave  length  ' . 
dep. 

Angle  dep. 
Angle  dep. 
Absolute  mag. 
Absolute  mag. 
Wave  length 
dep. 

Vifave  length 
dep. 

Absolute  mag. 

Angle  dep. 
Absolute  mag. 
Absolute  mag, 

V/ave  length 
dep. 

Angle  dep. 
Absolute  mag. 
Absolute  mag. 


Thickness  Ave.  Thlck- 
(cm.  X 10^)  ness (om.xlOl 

1.40±0.2  '‘7 
1.35±0.2  ( 

1.16^0.08  1.40 

1,25±0.07  r 

1.54±0.2  ) 


.74-0.15 

.70-0,4 

.90^0.07 

.82’^ 

.85-0.12 


1.00±Q.15 

1.32 

.90±0.2^' 

1.24±0.10 

1.18*0.08 


,9020.15 

1.04^0.08 

.67^0,07 

,80±0.08 


Ethylene  Glycol  16.7 


Absolute  mag,  >1*1 

Absolute  mag,  >2 
Absolute  mag,  >1.9 
Absolute  mag,  >1.6 


0.80 


1.16 


1.11 


0,85 


* Two  records 
**  One  record 
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mental  values.  This  scatter  was  relatively  large,  possibly  re- 
flecting uncertain  conditions  at  the  liquid-gas  interface  at  the  in- 
stant of  arrival  of  the  shock.  The  probable  errors  of  reflectivity 
values  ranged  from  11^  to  with  and  average  value  of  l5%.  The 
scatter  in  gas  calibration  reflectivities  yielded  a probable  error 
of  only  1,6^, 

TABLE  II 

It  will  be  noted  that  while  there  appears  to  be  some  discrep- 
ancy between  the  results  of  the  various  methods  of  measuring  L,  in 
most  cases  the  limits  of  error  overlap  for  identical  shocks.  The 
experimental  results  are  plotted  in  Figures  ij.,  5 and  6 together  with 
the  corresponding  theoretical  curves  computed  from  Lq,  (9), 
Conclusions 

The  comparison  between  the  experimental  and  theoretical  results 
for  acetone  and  ethyl  ether  suggests  that  either  the  shear  viscosity 
is  markedly  greater  under  these  conditions  or  that  a bulk  viscosity 
effect  is  present,  since  the  measured  values  of  L all  lie  signifi- 
cantly above  the  predicted  values.  The  values  of  nsed  in 

the  theoretical  computations  are  so  questionable,  however,  that  no 
clear  cut  conclusions  may  be  drawn  regarding  the  magnitude  of  the 
effect.  The  potentiality  of  the  method  for  gaining  information  as  to 
the  viscous  behavior  of  liquids  has,  however,  been  demonstrated;  an 
expression  for  shock  thickness  in  liquids,  valid  for  weak  shocks, 
has  boon  derived  from  hydrodynamic  theory,  and  an  apparatus  to 
measure  the  thickness  of  such  shock  waves  has  boeri  described.  The 
authors  wish  to  thank  Dr,  Edward  P.  Greene  for  valuable  help  and 
suggestions  throughout  the  experimental  work  and  Mr,  Thomas  L, 

Dawson  for  help  in  designing  and  constructing  the  apparatus. 
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